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ABSTRACT 

We investigate the behaviour of dissipative accreting matter close to a black hole 
as it provides the important observational features of galactic and extra-galactic black 
holes candidates. We find the complete set of global solutions in presence of viscosity 
and synchrotron cooling. We show that advective accretion flow can have standing 
shock wave and the dynamics of the shock is controlled by the dissipation parameters 
(both viscosity and cooling). We study the effective region of the parameter space for 
standing as well as oscillating shock. We find that shock front always moves towards the 
black hole as the dissipation parameters are increased. However, viscosity and cooling 
have opposite effects in deciding the solution topologies. We obtain two critical cooling 
parameters that separate the nature of accretion solution. 

Key words: accretion, accretion disc - black hole physics-shock waves. 



1 INTRODUCTION 

Accretion process around compact objects has been exten- 
sively studied during last three decades. Several attempts 
have been made for the development of the theory of ac- 
cretion disc. The standard theory of thin accretion disc 
provides a self-consistent s olution of Keplerian disc model 
flShakura fc Sunvaevlll973l . hereafter SS73). Soon after the 
standard disc model was proposed people realized that the 
disc are not Keplerian everywhe re, particularly in the radia- 
tion pressure domin ated region ijLiehtman fc Eardlevlll974l 
IChakrabart1ll996al) . In addition, this model (SS73) was un- 
successful to explain the origin of observed high energy ra- 
diatio n. A possible solution imme diately came in the liter- 
ature iSunvaev fc TitarchukllT985l) that disc may possess a 
Compton cloud which must inverse Comptonize the Keple- 
rian soft photons to produce the h ard X-rays. Incidentally , 
IChakrabarti fc Moltenil J1995T) and lLanzafame et alJ <ll998t) 
showed by numerical simulation that angular momentum 
distribution purely depends on the viscosity prescription un- 
der consideration and remains sub-Keplerian close to the 
black hole as flow must cross the horizon super-sonically 
JChakrabartil99d) . Recent observations do support the fact 
that sub-Keplerian matter may be pr esent in accretion disc 

JSmith et al.ll200lllSmith et alJl2002ft . 

In an earlier study, Chakrabarti llChakrabartil Il989t 
IChakrabartill99(Jl extensively discussed the transonic prop- 
erties of inviscid, polytropic flows and showed that standing 
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shock waves can form for a large region of parameter space 
spanned by the specific energy an d specific angular momen- 
tum of the flow. Vario u s authors llNoj>uta_fc_JIanawa| 



Yang fc Kafatol 11991 llu fc Yuan! Il99l iGu fc Luf 



Fukumura fc Tsurutal 120041) found the existence of stand- 
ing shock waves in different astrophysical contexts. In an 
accretion process around a black hole, flow is sub-sonic at a 
large distance and crosses the horizon super-sonically. Thus, 
black hole accretion is necessarily transonic. For a given set 
of input parameters, namely specific energy and specific an- 
gular momentum, the flow passes through the outer sonic 
point (xout) and remains super-sonic as it proceeds inwards. 
At a distance x s (< x ou t), flow is virtually stopped due to 
the resistance offered by the centrifugal barrier and conse- 
quently a shock is formed. The post-shock flow momentarily 
becomes sub-sonic and flow temperature increases in this 
region as the kinetic energy of the flow is converted in to 
thermal energy. Subsequently, flow picks up its velocity as 
it approaches black hole horizon and eventually crosses the 
horizon supersonically after passing through the inner sonic 
point (Xi n ). 

The first global and fully self-consistent solution of vis- 
cous transonic flows including advection in the optically thin 
and t hick limit was obtained by Chakrabarti dChakrabartil 
Il990l and references there in) where it was assumed that 
Keplerian flow at the outer edge can be come sub-Keplerian 
at the inner part of the disc. Later, IChakrabarti fc Dasl 
(2004a) presented a complete classification of global so- 
lutions of viscous transonic flow according to solution 
topologies and identified the region in the parameter space 
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where the flow possesses multiple sonic points. Within 
this region, if the viscosity parameter is below its criti- 
cal value (an < cms), there are solutions where Ra nkine- 
Hugoniot shock conditions jLandau fc Lifshit d Il95sl . here- 
after RHCs) are satisfied and standing shock is formed 
due to centrifugal barrier iChakrabartilll990l IChakrabartil 
Il996bl. IChakrabarti fc Dasl l2004al) . When viscosity is in- 
crease d beyond a critical value, shocks disappear (see Fig. 
12 of IChakrabarti fc DasH2004aD . Of course, the viscosity 
parameter strongly depends on the model parameters such 
as sonic point location and angular momentum at the inner 
boundary. T he importance of critical viscosity is being re- 
considered bv lGu fc Lul (l2004ll where a 'qualitative' estimate 
of critical viscosity for standing shock (an s ^ 0.5) was re- 
ported. Furthermore, many authors pointed out that shocks 
may underg o either radial or vert i cal or non-axisymmet rical 
oscillations jMolteni et alJ ll99Sl iGu fc Foriizzol l2003h al- 
though t hey are generally never r uined by these oscilla- 
tions. In iDas fc Chakrabartil ll2004h . Bremsstrahlung cool- 
ing also added keepi ng in mind that it is a very ineffi - 
cient cooling process llChattopadhvav fc Chakrabartil2000T) 
and complete set of global solutions of viscous transonic 
flow with and without shocks were presented. The hot and 
dense post shock flow which basically acts as a 'bound- 
ary layer' of the black holes could be the natural site of 
the hot radiation in the accretion disc and is considered to 
be a powerful tool in understanding the important astro- 
physic al phenomena such as the spectra l properties of black 
holes ( Chakrabarti fc Titarchukl Il995l IChakrabarti et alJ 
Il996d. Ebisawa et alJll996t iMandal fc Chakrabartil l200Eft . 
the source of quasi-periodic oscillation ( QPO) of the 
hard X-rays from the b lack hole candidates l|Molteni et alJ 
119961. iRvu et alJ Il997t IChakrabarti et alJ l2004bl) and the 
generation of accretion-pow ered relativistic bipolar out- 
flows /jets_j|D^^^L|j20^^_md references there in). Re- 
cently, IChakrabarti fc Mandall <l200 6) reported that this 
boundary layer could be responsible for the spectral state 
transition in at lea st some of the black h ole c andidates, such 
as Cy g XI. Also, iFender et aD (l2000l) and iDhawan etafl 
( 2000) showed that outflow could be produced from the same 
region that emits the Comptonized photons. A globally com- 
plete Infl ow-Outflow solution of viscous accretion flow is also 
found bv lChattopadhvav fc Dasl <2006h . 

Being motivated by the above arguments, in the present 
paper, we wish to follow up the stationary, axisymmetric, 
viscous accretion solutions around a Schwarzschild black 
hole in presence of Synchrotron cooling. Such dissipative vis- 
cous accreting flow have not yet been explored in the litera- 
ture. The space-time geometry around a non-rotating black 
hole is satisfactorily described by pseudo-Newtonian poten- 
tial dPaczvhski fc Wiitall980h We con sider similar viscosity 
prescription as in Tchakrabartil l)l996bF) . We calculate all rel- 
evant dynamical and thermodynamical flow variables and 
extensively study the dependence of such variables on the 
the flow parameters like specific energy, specific angular mo- 
mentum and the dissipation parameters (Synchrotron cool- 
ing and/or viscosity parameters). We identify the solution 
topologies which are essential for shock formation. In the 
present study, we concentrate on the non-dissipative shocks 
that preserve the energy across it. In a generalized accretion 
flow, cooling reduces the energy of the flow while viscosity 
not only tends to heat the flow but lowers the flow angular 



momentum as flow accretes. Thus the effect of viscosity and 
Synchrotron cooling on the dynamical structure of the accre- 
tion flow are expected to be different. Since shock has obser- 
vational consequences as mentioned above, it is, therefore, 
pertinent to know how the dynamics of shock is controlled 
by the dissipation parameters. In particular, whether shock 
front moves inward or outward if the dissipation is increased. 
In the present paper, we attempt to discuss this issue exten- 
sively. Moreover, we find that shock waves, standing or os- 
cillating, are still formed even when at very high dissipation 
limit. In this paper, we provide a useful formalism to study 
the dissipative transonic accretion flow which is valid for a 
wide range in accretion rates. 

We arrange the paper in the following way: in the next 
section, we present the governing equations. In §3, we per- 
form sonic point analysis. In §4, we study the global solu- 
tion topologies and classify the parameter space for solutions 
having multiple sonic points as a function of cooling param- 
eter. In §5, we study the dynamics of shock and dependence 
of shock properties on the the flow parameters. In §6, we 
classify the the parameter space as a function of cooling pa- 
rameter in terms of whether shocks, standing or oscillating, 
will form or not. In §7, we quantify the critical cooling pa- 
rameters. Finally, in §8, we make concluding remarks. 



2 BASIC HYDRODYNAMIC EQUATIONS 

We begin with a stationary, thin, viscous, axisymmetric ac- 
cretion flow on to a Schwarzschild black hole. In this pa- 
per, we consider ps e udo-N ewtonian potential introduced by 
IPaczvriski fc Wiital lll98CT) to solve the problem instead of 
using full general relativity, which allows us to use the New- 
tonian concept at the same time retaining all the salient 
features of the space-time geometry around a non-rotating 
black hole. The governing equations for the accreting flow 
are written on the equatorial plane of the accretion disc. The 
flow equations are made dimensionless by considering unit of 
length, time and the mass as r g = 2GMbh / '<? , 2GMbh / c? 
and Mbh respectively where, G is the gravitational con- 
stant, Mbh is the mass of the black hole and c is the ve- 
locity of light respectively. Henceforth, all the flow variables 
are expressed in geometrical units. 

In the steady state, the hydrodynamics of the axisym- 
metric accreting matter around a non-rotatin g black hole 
in the pseudo-Newtonian limit is given by (IChakrabartil 
Il996bft . 

(i) the radial Momentum equation : 



§ d& ldP _ ( . 

dx p dx x s ' 

(ii) the mass flux conservation equation : 

M = E*r 

(iii) the angular momentum conservation equation 

„d\(x) Id, 2 „, s. 

and finally, 

(iv) the entropy generation equation : 



ds 
dx 



hi) 



dP 

dx 



-yP dp 
p dx 



Q 



(i) 

(2) 
(3) 

(4) 
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The local variables x, p, P and A in the above 
equations are the radial distance, radial velocity, density, 
isotropic pressure and the angular momentum of the flow 
respectively. Here, 4>(x) = — denotes the pseudo- 
Newtonian potential and 7 is the adiabatic index of the flow. 
Furthermore, s and T are the specific entropy and the local 
temperature of the flow, E is the vertically integrated density 
and W x <t> is the viscous stress. Here, Q + and Q~ represent 
the energy gained and lost by the flow respectively, and M 
is the mass accretion rate. In our model, the accreting flow 
is assumed to be in hydrostatic equilibrium in the vertical 
direction. The half thickness of the disc is then obtain by 
equating the vertical component of the gravitational force 
and the pressure gradient force as, 



h(x) 



1/2 



(x-1). 



(5) 



where, a de notes the adiabatic sound speed defined as 
a ~ \flP~f Pi I n this paper, we us e the similar viscosity pre- 
scription of IChakrabartil lll996bl) . In this prescription, vis- 
cous stress remain continuous across the axisymmetric shock 
wave for flows with significant radial motion. 

Now we simplify the entropy equation (Eq. 4) as: 



ideal equation of state. In an electron-proton plasma, we use 
equipartition magnetic field (Eq. 7) to account for the non- 
dimensional synchrot ron cooling effect which is o btained in 
a convenient form as ( Sh apiro fc Teukolskvl l7983'): 



§x 3 / 2 {a 



with 



S = 7.5424 x 



Vf3Mp 2 



mi7 5 / 2 



2GM BH c 3 



(8) 



(9) 



where, e and m e represent the charge and mass of 
the electron respectively. We estimate elect ron tem - 
perature from the expression T e — \J m e /m p T p 
iChattopadhvav fe Chakrabartil [2002) and use it to 
obtain Eq. (8). Here, j3 is a dimensionless cooling parameter 
that controls the efficiency of cooling. When f3 — > 0, the 
flow is heating dominated as it cools most inefficiently 
In the present paper, we consider 7 = 4/3, r\ — 0.1 and 
M — 0.233MEdd for IOMq black hole, until otherwise 
stated. 



7- 



1 dP 

p dx 



■yP dp 
p 2 dx 



Q~ - 
ph 



A - T, (4a) 



where, r(= Q + / ph) denotes the energy gained due to the 
viscous heating and is obtained as, 



-ctnlnX / 2 . Q 2\ dfl 



(6) 



Here, an is the viscosity parameter and Q.(x) is the 
angular velocity of the accreting matter at a radial dis- 
tance x. The polytropic index is denoted by n = (1 — 7)" 1 - 
I n and I n +i come from the vertica lly averaged density 
and pressure iMatsumoto et ail I1984T) and g — I n +i/I n - 
The term A(= Q~/ph) represents the energy loss by the 
flow. In the present analysis, we ignore energy loss due to 
Bremsst rahlung cooling as it is not very effi cient cooling 
process dChattopadhvav fc Chakrabartil I2OOOI) and include 
synchrotron cooling only. Indeed, magnetic field is ubiqui- 
tous inside the accretion disc and therefore, the ionized flow 
in turn should emit synchrotron radiation that may cool 
down the accreting flow significantly. Since the satisfactory 
description of magnetic field inside the accretion disc is not 
well understood, we thus assume only random or stochastic 
magnetic field. Such magnetic field may or may not be in 
equipartition with the accretion flow. We define a parame- 
ter r\ as the ratio of the magnetic pressure and the thermal 
pressure of the flow under consideration and is given by, 

V = ff^r (7) 

where, B denotes the magnetic field strength, ks is the 
Boltzmann constant, p is the mean molecular weight and 
m v is the mass of the proton respectively. In general, 77 <: 1 
which ensures that magnetic fields definitely rema in con- 
fined within the disc jMandal fc Chakrabartal2005l) . While 
obtaining Eq. (7), we consider the gas is assumed to obey an 



3 SONIC POINT ANALYSIS 

We solve Eqs.(l-3, 4 a) following the s tandard method of 
sonic point analysis llChakrabartil ll989Tl . We calculate the 
radial velocity gradient as: 

di3 N 

Tx = D> ^ 
where, the Numerator iV is given by, 



N 



x 3 / 2 (x - 1) 



Qn-/n(a 2 g + 7^ ) a^gl ^j^x - 3)(o 



J 2 X 



A 



x 3 2(x - l) 2 



7 2 x(a; — 1) 
( 7 + l)tf 2 2a n p/ n (a 2 fl + 7 tf 2 )" 



(7-1) 



■ff 2 a 2 (5x - 3) 2\a n I n tf(a 2 g + -yd 2 ) 



£■(7 — l)(x — 1) 
and the denominator D is, 



73; z 



(10a) 



D 



2a 2 -Q (7+l)i9 3 a 2 J I„^(a 2 g + -yd 2 



(7 ~ 1) (7 - 1) 7 
The gradient of sound speed is obtained as 



(2ff-l)-^ 



a~ 
~$ 2 
(106) 



da fa ji3 \ di3 ~ 



dx 



■d a J dx a 



A 



2(3-1) 



+ 



(5a; — 3)a 



2x(x - 1) 
(11) 

and similarly, the gradient of angular momentum is calcu- 
lated as: 



d\ ctn(a 2 g + 71? 2 ) anx(a 2 g — 'yd 2 ) d-d 2otjiaxg da 

dx 71? 71? 2 dx 71? dx 

(12) 

At the outer edge of the accretion disc, matter starts 
to accrete with almost negligible radial velocity and enters 
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into the black hole with velocity of light. This indicates that 
'radial velocity gradient' should be always real and finite to 
maintain the accretion flow smooth along the streamline. 
However, Eq. (10b) indicates that there may be some points 
where denominator (D) vanishes. Since the flow is smooth 
everywhere along the streamline, the point where denom- 
inator tends to zero, the numerator (JV) must also vanish 
there. The point where both the numerator and denomina- 
tor vanish simultaneously is called as critical point or sonic 
point. Setting D = 0, one can easily obtain the expression 
for the Mach number (M — •d/a) at the sonic point as, 



/ -mo — \ m« — 4mim3 
M(x c ) = \ V 2 (13) 



2mi 



where, 



mi = 7 2 [a 2 n I n (2g - 1)( 7 - 1) - ( 7 + 1)] 
m 2 = 2 7 [7 + anl n g{g - 1)( 7 - 1)] 

m 3 = -anl n g 2 (l - 1) 

In the weak viscosity limit, the Mach number at the 
sonic point becomes, 



M(x c ) = 



7 + 1 



(13a) 



This result is exactly same as IChakrabartil il989T) . This in- 
dicates that Mach number at the sonic point in presence of 
different cooling process remains same as in the case for the 
non-dissipative accretion flow. We obtain an algebraic equa- 
tion for sound speed by using another sonic point condition 
iV = which is given by, 

F(£ c ,X Cl x c ) = Aa 3 {x)+Ba(x)+Ca{x)+V = 0, (14) 
where, 

S 



A= — 



x 3 / 2 (x- 1)' 



(14a) 



B = ahl n (g + jM 2 ) 2 a^I n g{hx - 3)(<? + jM 2 ) 



"f 2 X 



y 2 x(x — 1) 



M 2 (5a:-3) 
x (l — l)( x ~ 1) ' 

2\a n InM(g + 7M 2 ) 
7X 2 



(146) 
(14c) 



and 



A 



x 3 2{x - l) 2 



( 7 + l)M 2 2a 2 lg J n (ff + 7 M 2 ) 
(7-I) 7 



W e solve Eq. (14) analytically iAbramowitz fc Steeunl 
Il97(t) and obtain the sound speed at the sonic point by 
knowing the flow parameters. Depending on the input pa- 
rameters, a non-dissipative flow may have maximum of four 



sonic points iDas et al.l2001bT) . Among them, one always lies 
inside the horizon and the others, if they at all exist, remain 
outside the horizon. In our present analysis, we continue to 
use similar concept to obta in dissipative accretion solutions 
JChakrabarti fc Dasl feotMa). The nature of sonic points and 
its detailed properties could be easily understood with the 
extensive study of Eq. (10) when the flow parameters are 
known. Infact, nature of sonic points depends on the value 
of velocity gradients at the sonic point. In general, d-d/dx 
possesses two values at the sonic point. One of them is for 
accretion flow and the other is valid for wind. When both the 
derivatives are real and is of opposite signs, the sonic point 
is saddle type or X-type. Nodal type sonic point is obtained 
when both the derivatives are real and is of same sign. If 
the derivatives are complex, the sonic point belongs to the 
class of spiral type or O-type. In the astrophysical context, 
saddle type sonic point has some special importance since 
transonic flows usually pass through it. In addition, in or- 
der to form a standing shock, the flow must have more than 
one saddle type sonic points among them the closest and 
furthest values corresponds to inner and outer X-type sonic 
points. The O-type sonic point which lies between the two 
saddle type sonic point is unphysical in the sense that no 
stationary solution passes through it. 



4 SOLUTION TOPOLOGIES 

To obtain a complete solution, one needs to supply the 
boundary values of energy, accretion rate and angular mo- 
mentum of the flow for a given viscosity (an) and cooling 
parameter (/3) and integrates Eqs.(10, 11, 12) simultane- 
ously with the help of Eqs. (10a, 10b). Since black hole so- 
lution is necessarily transonic, flow must pass through the 
sonic points. In addition, flow presumably originates from 
the Keplerian disc (xKep) and it must carry angular mo- 
mentum X(xKe P ), a part of which is transported due to vis- 
cosity outwards and the other part is advected into the black 
hole through the sonic point. The distribution of the angular 
momentum from the inner edge of the disc to the XKep is 
obtained from Eq. (12). Incidentally, the acceptable range 
of angular momentum at the inner edge of the disc and the 
inner sonic point location of the accreting sol ution is very 
small (say, 1.5 < Xin < 2 and 2 < Xi„ < 4, IChakrabartil 
Il989t IChakrabarti fc Dasl l2004al) . 'We. thus, choose inner 
sonic point location Xi n and the angular momentum at the 
inner sonic point X(xi„) (hereafter, we denote it as Xin) as 
input parameters instead of the aforementioned energy and 
the angular momentum flow at the outer edge of the disc. We 
once integrate Eqs. (10, 11, 12) inwards upto the horizon and 
again outwards to a large distance (~ XKep) to obtain a com- 
plete transonic solution. A comprehens ive study of viscous 
transon ic flow is already presented by IChakrabarti fc Dasl 
(2004a) and therefore, in the current paper, we mainly con- 
centrate on the properties of accreting solution in presence 
of synchrotron cooling, until otherwise stated. 



4.1 Behavior of Global Solutions 

In order to obtain a shock solution around a black hole, ac- 
creting flow must possess multiple saddle type sonic points. 
Shock joins different solution branches, one passes through 
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Figure 1. Global solution topologies in presence of Synchrotron 
cooling. In each boxes, variation of Mach number is plotted with 
logarithmic radial distance. Flow parameters are Xi n = 3.107, 
^in = 1-6. Cooling parameters are marked in the figure. 

the inner sonic point and the other passes through the outer 
sonic point. Thus, it is important to understand the nature 
of solution trajectories in presence of energy dissipation. 

Fig. 1 shows examples of global solutions of inviscid ac- 
cretion flow passing through the inner sonic point. Each box 
depicts Mach number (M = l?/a) variation as a function of 
the logarithmic radial distance. We write the cooling param- 
eter in each box (marked as 0.0, 0.05, 0.1 and 0.2). The 
other flow parameters for Fig. 1 are chosen as x c — 3.107, 
Xi„ = 1.6 and n = 3 respectively. Here, all the sonic points 
are saddle type, and the solutions may or may not have shock 
depending on its global solution topologies. In the first col- 
umn, the flow is dissipation free and the sonic points are 
uniquely determined as energy and angular momentum are 
conserved all throughout. Accreting solutions of this kind 
deviate from the Keplerian disc at the outer part of the disc 
and enter into black hole super-sonically with the radial ve- 
locity equal to the velocity of light after crossing the outer 
sonic point. Another possibility is that super-sonic flow may 
suffer shock transition at sub-sonic branch, if shock condi- 
tions hold (see §5) and crosses the black hole horizon after 
passing through the inner sonic point as it has to satisfy the 
'super-sonic' inner boundary condition at the horizon. In the 
second column, cooling is incorporated and the flow topolo- 
gies passing through the inner sonic point start to open up. 
The outer sonic point, if it exits, can not be determined un- 
less a shock is formed as energy of the flow is not conserved 
here. When cooling is increased further, the flow topologies 
are completely opened up (third and fourth column) and 
they leave the accretion shock regime. In this particular sit- 
uation, flow deviates from the Keplerian disc far away from 
the black hole and enters into it straight away through the 
inner sonic point. So, for a given set of flow parameters, 
there must be a critical cooling parameter (/3 CC ) for which 
closed topologies passing through the same inner sonic point 
become open topologies. An exhaustive discussion of critical 
cooling parameter will be presented in §7. 

4.2 Description of the parameter space 

In order to understand the dissipative accretion flow, here, 
we classify the parameter space as a function of the cooling 
parameter in the £i n -A in plane where £j n denotes the en- 
ergy of the flow at the inner sonic point (xj n ). Infact, in a 
cooling dominated inviscid accretion flow, angular momen- 
tum remain conserved along the streamline of the flow. In 
Fig. 2, we have identified the region of the parameter space 
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Figure 2. Classification of parameter space for closed topolo- 
gies passing through the inner sonic point. The region bounded 
by different curve is obtained by using different cooling factor 
[P = 0.0 (solid), = 0.00787 (dotted), f3 = 0.0787 (dashed) and 
P = 0.393 (dot-dashed)]. Notice that, parameter space gradually 
shrinks with the increase of cooling factor. 

for closed spiraling accretion solution which passes through 
the inner sonic point as in panel 1-2 of Fig. 1. Accretion 
flows of this kind possess multiple sonic points and flow 
may suffer shock transition if the outer sonic point exits. 
The shock would be stationary or oscillating depending on 
whether the Rankine-Hugoniot relation is satisfied or not. 
The region bounded by the solid curve is obtained for dis- 
sipation free (0 = 0) acc r etion flow and it coincides with 
the result of IChakrabartil lll989t) . The dotted, dashed and 
the dot-dashed regions are obtained for cooling parameters 
= 0.00787, 0.0787 and 0.393 respectively. For higher 0, the 
parameter space for closed topologies passing through the 
inner sonic point reduces in both lower and higher angular 
momentum sides. This is not surprising, since an increase of 
cooling induces damping inside the flow which reduces the 
number the sonic points for the same set of initial parameter 
as seen in Fig. 1. With the rise of 0, the effective region of 
parameter space for the closed topologies gradually shrinks 
and disappears completely when the critical cooling (0 CC ) is 
reached. 



5 SHOCK SOLUTIONS AND 

CLASSIFICATION OF PARAMETER SPACE 

The existence of shock waves in accretion flow has been re- 
ported in many astrophysical context and the thermody- 
namical properties of shock waves like its location, strength 
and compression ratio could be determined exactly by using 
the RHCs. In the present stud y, we continu e to use sim- 
ilar viscosity prescription as in IChakrabartil lll996bl) that 
maintains viscous stress {W x $) to be continuous across the 
shock in presence of significant radial motion. We compute 
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Figure 3. A complete accretion solution topology drawn with 
standing shocks (x s = 23.83). See text for details. 

the shock locations following an usef ul technique first pre- 
sented bv IChakrabarti fe Dasl i2004al) for dissipative accre- 
tion flows. 



5.1 Shock Dynamics and Shock properties 

In Figure 3, we present a complete solution of vertically 
averaged flow where a shock wave connects two solution 
branches — one passing through the outer sonic point (O) 
and the other passing through the inner sonic point (I). In 
this figure, we plot the Mach number variation with the 
logarithmic radial distance. The flow parameters are cho- 
sen as x in = 2.5308, A in = 1.725 and j3 = 0.0157 respec- 
tively. Sub-sonic accretion flow crosses the outer sonic point 
at x out = 80.07 and makes discontinuous transition (shock) 
at x s = 23.82 before entering into the black hole. The shock 
location is shown by vertical dotted line. The exact location 
of the shock is obtained by solving Rankine-Hugotinot re- 
lations. In the present study, shock is assumed to be thin 
and non-dissipative. In addition, we ignore the presence of 
any excess sources of torque at the shock, that maintains the 
continuity of angular momentum across it. At the shock, en- 
tropy is generated which is subsequently advected into the 
black hole. 

In Fig. 4, we have shown how the shock location (indi- 
cated by the vertical lines) changes with the dissipation pa- 
rameters (an and/or f3) for a given set of initial flow param- 
eters. We inject matter sub-sonically at the outer edge of the 
disc Xi„j = 145 with local energy £(x in j) = 3.3663 x 10~ 3 
and angular momentum X(xi„j) = 1.725. (hereafter, we de- 
note them as £i n j and Xinj respectively). The sub-sonic flow 
first passes through the outer sonic point and becomes super- 
sonic. The Rankine-Hugoniot condition are satisfied here 
and eventually a shock is formed. The solid vertical line rep- 
resents the shock location (x s = 25.99) for non-dissipation 
(an = and j3 = 0) flow. As viscosity is turn on, shock front 




log(x) 

Figure 4. Plot of Mach number with the logarithmic radial dis- 
tance. Flows with same energy £i„j = 3.3663 X 10 -3 and angular 
momentum Xi n j = 1.725 at the outer edge Xi n j = 145 are injected 
with different dissipation parameters. Dissipation parameters are 
marked in the figure. The corresponding shock locations are indi- 
cated by the vertical lines [x s = 25.99 (solid), x s = 20.6 (dashed), 
x a = 22.94 (dotted) and x a = 18.62 (dot-dashed),]. 

(a; s = 20.6) moves forward denoted by the dashed vertical 
line. Due to viscosity, accreting flow losses angular momen- 
tum causing the reduction of centrifugal pressure and at the 
same time energy of the flow increases. At the shock, total 
pressure must be continuous and thus, the dynamics of the 
shock are mainly controlled by the resultant pressure across 
it. Since the shock moves in and as X(x) get reduced along 
the viscous flow, therefore, we can conclude that centrifugal 
force is the primary cause for shock formation. When syn- 
chrotron cooling is effective in an inviscid accretion flow, 
shock location (x s — 22.94) again proceeds towards the 
horizon depicted by vertical dotted line. In the pre-shock 
region, density and temperature are low and therefore, the 
effect of cooling is negligible (see below). The situation is 
completely opposite in the post-shock region and cooling 
reduces the post-shock pressure causing the shock to move 
forward further. This indicates that thermal pressure is also 
not less important in determining the final shock location. In 
a generalized accretion flow where both viscosity and cool- 
ing are present, shock location is predicted at x a = 18.62 
for the same set of input parameter as above and is de- 
noted by vertical dot-dashed line. In this particular case, 
shock front is shifted significantly due to combined effects 
of viscosit y and cooling. When the ma ss loss from the vis- 
cous disc Ichattopadhvav fc DaJ F2006) is considered, shock 
is seen to form even closer to the black hole as some part of 
the accreted matter is ejected as outflows which reduces the 
post-shock pressure drastically. 

In Fig. 5a, we have plotted the logarithmic proton tem- 
perature (T = /j,m p a 2 I^b) of a weakly viscous flow as a 
function of the radial distance for the same set of input 
parameters as in Fig. 4 with (3 = 0.0236. Clearly Fig. 5a 
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Figure 5. The variation of (a) Flow temperature and (b) total 
Synchrotron loss (in geometric units) with radial distance. 



indicates that post-shock temperature is much higher com- 
pared to the pre-shock temperature. This is because shock 
compresses the flow and makes it hotter. In particular, at 
the shock, flow temperature rises sharply as the kinetic en- 
ergy of the flow converted to the thermal energy there. In 
Fig. 5b, we show how the total synchrotron loss (in logarith- 
mic scale with geometric unit) vary with the radial distance 
for the same set of flow parameters as in Fig 5a. The net 
energy loss is negligible in the pre-shock region while it in- 
creases significantly as the post-shock flow advected towards 
the black hole. 

Fig. 6a compares the shock locations as a function of the 
cooling parameter /3 for different values of angular momen- 
tum Xinj- In all the cases, the matter is injected at the outer 
edge of the disc x in j = 145. The solid curve is for \ in j = 1.7 
and the dotted and dashed curve are for Xi„j = 1.725 and 
Xinj = 1-75 respectively. The corresponding local energies 



at the injection point are £i„ 



2.5156 x 10- 3 (solid), 



3.3663 x 1CT 3 (dotted) and 4.3474 x 1CT 3 (dashed) respec- 
tively. This example shows that stable shocks can form for 
a wide range of cooling parameter. For a given local energy 
and angular momentum of the flow at the injection point, 
shock location decreases with the increase of the cooling pa- 
rameter as flow losses energy as it accretes. When cooling 
is above its critical value (/3 > /3 CS ), shock ceases to exist 
as the RHCs are not satisfied there. An important point to 
note is that the shock forms at larger radii for higher angu- 
lar momentum as the centrifugal barrier becomes stronger 
with angular momentum. 

It is useful to study the density distribution across 
the shock since cooling processes as well as emergent 
radiation from the disc ar e strongly depends on it 
llChakrabarti fc Mandall l2006). For this we compute com- 
pression ratio (R) which is defined as the ratio of the ver- 
tical averaged post-shock density to the pre-shock density. 
As an example, in Fig. 6b, we show the variation of com- 
pression ratio with the cooling parameter for the same set 
of flow parameters as in Fig 6a. For a given set of local en- 
ergy and the angular momentum of the flow, compression 
ratio R increases monotonically with the cooling parame- 
ter. In a convergent accretion flow, shock is pushed inwards 
when cooling is enhanced. This causes more compression in 
the post-shock region. There is a cut-off at a critical cooling 
limit as shock disappears there. 

In our further study of shock properties, in Fig. 6c, we 
draw the variation of the shock strength (O) with the cool- 
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Figure 6. (a) Variation of shock location with the cooling fac- 
tor when the flows with identical angular momentum and en- 
ergy are injected from a fixed outer boundary. Solid, dotted 
and dashed curves are drawn for flow with angular momentum 
^in = 1.7, 1.725 and 1.75. Shock forms further for higher an- 
gular momentum flow suggests that shock are mainly centrifu- 
gally driven. In all cases, shock moves closer to the black hole as 
the cooling factor is increased, (b) Variation of compression ratio 
(R = S+/E_) with the cooling factor for same set of parameters 
as (a). Subscripts "+" and "-" refer to quantities before and after 
shock. Compression ratio monotonically increases in all the cases 
as the cooling factor is increased, (c) Variation of shock strength 
(0 = M—/M+) with the cooling factor for same set of parameters 
as (a). Shock strength gradually increases in all the cases as the 
cooling factor is increased. 



ing parameter. Shock strength is defined as the ratio of the 
pre-shock Mach number to the post-shock Mach number and 
it determines the temperature jump at the shock. Same set 
of initial parameters as in Fig. 6a are used here. In the dis- 
sipation free limit, the strength of the shock is weak and it 
becomes stronger as cooling parameter is increased. 

5.2 General behavior of Shocks 

So far, we only study the properties of shock waves when 
matter is injected from a fixed outer edge of the disc. In 
reality, global solutions that contain shock waves are not 
isolated solutions, but are present in a large range of the 
flow parameters. We continue our study of shock properties 
for a wide range of flow parameters. In Fig. 7a, we draw the 
shock location (x s ) as a function of specific energy at the 
inner sonic point The different curves are drawn for a 
set of cooling parameters starting form (5 = 0.0 (rightmost) 
to 0.144 (leftmost) with an interval A/J = 0.016. Here, the 
angular momentum at the inner sonic points is chosen as 
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Figure 7. (a) Variation of shock location with flow energy at 
the inner sonic point. The different curves are drawn for different 
cooling factor. The rightmost curve is drawn for dissipation free 
flow. As dissipation increases (A/3 = 0.016), shock forms at a 
closer distance, (b) Variation of compression ratio as a function 
of flow energy at the inner sonic point and the cooling factor. The 
flow parameters are same as (a). Compression ratio decreases in 
both higher and lower ends as the cooling is enhanced. 



Figure 8. Modification of Parameter space for standing shock 
as a function of cooling parameter. Parameter space shrinks with 
the increase of cooling factor. 

viously expected to be produced in presence of synchrotron 
cooling. 



\i n = 1.75. Notice that shock forms only for a particular 
range of £i n and in the weak cooling limit, shock location is 
varied in a wide range (from 10 - 80 r g in this particular set 
of parameters). At higher cooling parameter /3, shocks form 
at closer radii and the range of shock locations gradually 
shrinks. In particular, the upper limit of the shock locations 
proceeds more rapidly towards the black hole horizon with 
the increase of cooling parameter. When /3 > /3 CS is reached, 
accreting flow fails to satisfy the RHC and shock disappears. 
Shocks can even form at negative energies at the inner sonic 
point (in general, the flow is supposed to have positive en- 
ergy at the shock) indicates that due to cooling, flow looses 
significant amount of energy while falling into the black hole. 

In Fig. 7b, we plot the variation of the compression ra- 
tio with £i n for the same set of initial flow parameters used 
in Fig. 7a. Each individual curve is obtain for a given cooling 
parameter ranging from /3 = 0.0 (rightmost) to /3 = 0.144 
at intervals of A/3 = 0.016. Strongest shocks (R — * 7) form 
only in the weak cooling limit. Figure clearly shows that the 
upper and lower limit of compression ratio monotonically 
decrease for the gradual increase of the cooling parameter 
and finally compression merges t o the intermedia te value 
(R ~ 4, a measure of strong shock. iDas et al.ll2001al and ref- 
erences there in) upto a critical cooling limit beyond which 
shock does not form. Thus, strong shocks are formed for a 
large range of cooling parameter and outflow and jet are ob- 



6 PARAMETER SPACE FOR SHOCK 

In Fig. 8, we redraw the parameter space as in Fig. 2, but 
consider the formation of standing shock only in a cool- 
ing dominated flow. The parameter space is expected to be 
modified with the increase of cooling parameter. We iden- 
tify the modified regions of the parameter space for sta- 
tionary shock using various cooling parameter. We write 
cooling parameters in the figure. For /3 = 0.0, the region 
bounded by the so lid curve identically merges with that of 
IChakrabartil (^989). At higher cooling parameter, the effec- 
tive region of the parameter space for standing shock is re- 
duced in both the higher and lower angular momentum side 
and it is shifted to negative energy region as well. Beyond a 
critical cooling limit this region disappears completely. 

We continue our study of parameter space in the energy- 
angular momentum (£i„,\i n ) plane according to the accre- 
tion flow topologies for f3 = 0.00787 in the weak viscosity 
limit. The boundary denoted by the solid curve in Fig. 9 
identifies the region of the parameter space for closed topolo- 
gies passing through the inner sonic point (Fig. 1). Fur- 
ther sub-divisions of parameter space are marked by dashed 
and dotted curves depending on the behavior of the solu- 
tion topologies. Accretion topologies with parameters chosen 
from different region of the parameter space are presented 
in the small boxes at the bottom left of Fig. 9. In each small 
boxes, we plot the variation of Mach number as a function 
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Figure 9. Division of parameter space in the Ei n , Xi„ plane ac- 
cording to the topologies. Accretion solutions with parameters 
taken from different parts of the parameter space are shown in 
the insets. 
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Figure 10. Region of the parameter space for standing shock 
for different dissipation parameters. Dissipation parameters are 
marked in the figure. 



of logarithmic radial distance. Flows with parameters from 
the region ABD possess multiple sonic points. In particular, 
when the flow parameters fall in the region ABC separated 
by the dashed boundary, RHCs are satisfied. Accretion so- 
lution of this kind is drawn in box labeled S. But, flows with 
parameters from rest of the region of ABT), RHCs are not 
satisfied, however, the entropy of the flow at the inner sonic 
point continues to be higher compared to that at the outer 
sonic point. In this case, shock starts os cillating causing a 
periodic breathing of post-shock region IIRvu et alJ ll997). 
The box marked OS shows an accretion solution having os- 
cillating shock. For higher cooling, we o btain a new solution 
topolo gy as shown in box CIM as seen in lDas fc Chakrabartil 
(2004). We draw this solution with dotted curve as it is ob- 
tained for higher cooling parameter. This multi-valued solu- 
tion is expected to be unstable and may cause a non-steady 
accretion. The accretion solution with parameters from the 
I region of the parameter space is drawn in the box labeled 
I. This accretion solution straight away passes through the 
inner sonic point before entering into the black hole. So- 
lution inside the box O represents an accretion flow which 
passes through the outer sonic point only. The box marked 
CI shows an closed accretion solution passing through the 
inner sonic point. This type of solution does not extend to 
the outer edge to join smoothly with any Keplerian disc and 
flow is expected to be unstable. When dissipation is signif- 
icantly increased, the closed topology of CI opens up (Fig. 
1) and then flow can join with the Keplerian disc. 

Fig. 10 shows the comparative study of the parame- 
ter space for standing shocks in different dissipation limits. 
The dissipation parameters are marked. The solid bound- 
ary separates the parameter space for standing shock in 



the non-dissipative accretion flow llChakrabartill989Fl . When 
the viscosity is included, the effective region of the param- 
eter space for standing shock reduces and shifts towards 
the higher energy and the l ower angular momentum regime 
JChakrabarti fc Das! l2004al) . The region bounded by dot- 
dashed curve is obtained for standing shocks in viscous flow. 
On the other hand, for a cooling dominated inviscid ac- 
cretion flow parameter space for standing shock shrinks in 
both ends of the angular momentum (lower and higher) and 
moved to the lower energy sides (present paper). This pa- 
rameter space is drawn with dashed curve. When viscosity 
and cooling are both present, the parameter space for shock 
settles down at a intermediate region (shaded part). This 
indicates that viscosity and Synchrotron cooling act oppo- 
sitely for deciding the shock parameter space. However, they 
behave similarly while controlling the dynamics of the shock 
waves (Fig. 4). 



7 CRITICAL COOLING PARAMETER 

In our earlier discussion, we already mentioned that the 
global behavior of the cooling dominated accretion flow 
changes when the cooling parameter exceeds its critical 
value. In particular, there exists two critical cooling parame- 
ters. The first one identifies a boundary which separates the 
closed topologies from the open topologies passing through 
the inner sonic point. The second one is obtained at the re- 
gion of the closed solutions based on the criteria of whether 
standing shock is formed or not. These critical cooling pa- 
rameters of course strongly depend on the inflow parameters. 
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Figure 11. Variation of critical cooling parameter as a function 
of angular momentum at the inner sonic point. The region out- 
side the dotted curve allows solutions with open topologies and 
the region shaded by solid lines separating the standing shock 
solutions from the oscillating shock. 

In Fig. 11, we present the plot of critical cooling parameters 
against angular momentum at the inner sonic point. Differ- 
ent regions are marked. Figure shows that closed topologies 
as well as shocks are formed in the lower and higher angular 
momentum side for smaller cooling parameters. As cooling 
parameter is increased, shock forms only somewhere at the 
intermediate angular momentum domain which is consistent 
with our earlier discussion. 



8 CONCLUDING REMARKS 

In this paper, we have presented the properties of the viscous 
accretion flows around a stellar-mass black hole in presence 
of synchrotron cooling. In a realistic accretion flows, the en- 
ergy loss due to radiation is not negligible and it has sig- 
nificant effect on the dynamical structure of the flow as we 
obtain here. This will alter the emitted spectrum and lumi- 
nosity as well. 

In the past, efforts were made to study the properties of 
viscous transonic flow in the cool ing free limit dChakrabartil 
Il99d. IChakrabarti fc Das! l2004atl and /or cooling is consid- 
ered in an approximate way Ichakrabartll Il996fj) . No at- 
tempts were made to study a complete accretion solution 
which may contain shocks in presence of synchrotron cool- 
ing. Here, we have obtained a complete set of global solu- 
tions. We showed that for a large region of the parameter 
space a transonic flow can have shock waves when viscos- 
ity as well as synchrotron cooling is very high. We find that 



the shock locations may vary from around ten Schwarzschild 
radii to several tens of Schwarzschild radii depending on the 
inflow parameters. We showed that standing shocks where 
significant amount of kinetic energy to thermal energy con- 
version takes place can form much closer to the black hole 
with the enhancement of dissipation. For higher viscosity, 
flow transports angular momentum more and more, as the 
flow moves inward and at the same time it tends to heat the 
flow. This causes the centrifugal pressure to be much lower at 
the post-shock region although the post-shock thermal pres- 
sure is higher. When the viscosity parameter is increased, 
shock forms closer to the black hole horizon, indicates that 
perhaps shocks are centrifugal pressure supported. On the 
other hand, Synchrotron photons generated by the stochas- 
tic/random magnetic field can cool down the flow and in 
addition, the flow angular momentum distribution remains 
insensitive to it (Eq. 3). As the synchrotron cooling is in- 
creased, the post-shock flow is cooled down more efficiently 
and the thermal pressure in this region reduces significantly. 
This causes the shock to move inwards to maintain pressure 
balance (thermal plus ram) across it. Thus, viscosity and 
cooling play similar role as far as the dynamic properties 
of the shock is concerned. However, cooling can not com- 
pletely nullify the effect of viscous heating as they depend 
differently on the flow variables. 

We obtain the parameter space for standing shock waves 
as a function of cooling parameter and find that the effec- 
tive region of the parameter space shrinks as the cooling 
parameter is enhanced. We also identify the region of the 
parameter space where standing shock conditions (RHCs) 
are not favourable but accreting solutions having two saddle 
type sonic point are pres ent. These solut ions generally ex- 
hibits non-steady shocks dRvu et alll997tl and it could even 
form when cooling is very high. Moreover, the possibility of 
shock formations decreases for increasing cooling parameter 
and when cooling parameter is increased beyond a critical 
value, shocks disappear (Fig. 11). We also find that viscosity 
and Synchrotron cooling induce completely opposite effect 
in deciding the parameter space for stationary shock waves. 
Indeed, it is clear that cooling does affect the dynamics of 
the accretion flow as well as shock parameter s pace quite 
significantly. Therefore, cooling free assumption JGu fc Lul 
2004) is perhaps too simplistic in a general transonic accre- 
tion flow. 

We have discussed that the black hole accretion models 
may need to include shock waves as they provide a com- 
plete explanation of the o bserved steady state, as well as 
time dependent behaviour JChakrabarti fc Titarchuklll995l. 
iRvu et alJll997l) . In particular, as cooling is increased, post- 
shock energy goes down and shock location is reduced . 
In addition, earl ier studies llChakrabarti fc Manickamll2000l . 
IRvu et ail Il997f) show that quasi-periodic oscillations of 
hard radiations from black holes is proportional to the in- 
fall time and therefore, QPO frequency increases with the 
enhancement of cooling parameter. Observational findings 
do support this view that for higher accretion rate, QPO 
frequency increases. Therefore, possibly shocks may be an 
active ingredient in the accretion flow. 

In the present paper, we have shown that for a given 
set of flow parameters there exist two critical cooling pa- 
rameters which separates the parameter space into three re- 
gions (Fig. 11). In the first region, flow passes through a 
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single X-type sonic point. In region two, flow has multiple 
sonic points but Rankine-Hugoniot relations are not favor- 
able anywhere while in the third region Rankine-Hugoniot 
relations are satisfied and standing shocks form. 

An important point we need to discuss is that since the 
basic physics of advective flows are similar around a neu- 
tron star (differs from black hole solution only due to inner 
boundary conditions), the conclusion may remain roughly 
similar although the boundary of the neutron star, where 
half of the binding energy could be released, may be more 
luminous than that of a black hole. 

In our calculation, we have made some approximations, 
such as we consider stochastic magnetic field instead of large 
scale field. We ignore the effect of radiation pressure. Also, 
we assume a fixed polytropic index rather than computing it 
self-consistently. We use Paczyhski- Wiita pseudo-potential 
iPaczvriski fe W iita 1980) to describe the space-time geom- 
etry around a Schwarzschild black hole. However, we believe 
that our basic conclusion will not be affected qualitatively by 
these approximations. More generalized calculations around 
a rotating black hole is beyond the scope of the present pa- 
per and will be reported elsewhere. 
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